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We analyse pp elastic scattering data at the highest ac- 
celerator energy region (10 < -^s < 62.5 GeV) through a 
multiple diffraction approach. The use of Martin's formula in 
a model developed earlier is substituted by the introduction of 
a complex elementary (parton-parton) amplitude. With this 
the total cross section and the p-parameter may be simulta- 
neously investigated and, with the exception of the diffraction 
minimum at some ISR energies, a satisfactory description of 
all experimental data is obtained. Total cross sections ex- 
trapolations to cosmic-ray energies (-^s > 6 TeV) show agree- 
ment with the reanalysis of the Akeno data performed by N. 
N. Nikolaev and also with Gaisser, Sukhatme and Yodh re- 
sults, leading to the prediction af^^{^ = 16 TeV) — 147 mb. 
Physical interpretations and critical remarks concerning our 
parametrizations and results are also presented and discussed. 

13.85.Dz, 11.80.La, 11.80.Fv 



I. INTRODUCTION 

Elastic proton-proton scattering is the most simple 
process in high-energy hadronic interactions. Despite 
the amount of experimental data available and model 
descriptions of these data, a treatment based exclusively 
in the field theory of strong interactions, the quantum- 
chromodynamics (QCD), is still missing. Beyond the in- 
trinsic interest coming from this fact, new expectations 
are associated with the next proton collider, the CERN 
Large Hadron Collider (LHC), optimistically planned to 
begin experiments at 14 TeV in 2005. The situation sug- 
gests that a fundamental and difficult task for present 
and future developments, is to find connections between 
model descriptions and reliable calculational schemes in 
QCD. Since elastic scattering incorporates soft processes 
one expects that nonperturbative formalisms will play a 
fundamental role. 

Based on the above considerations we have investigated 
elastic proton-proton scattering through phenomenolog- 
ical approaches and, simultaneously, looking for connec- 
tions with nonperturbative QCD treatments. With this 
strategy, we improved some aspects of a multiple diffrac- 
tion model developed earlier. The main points concern 
the use of a complex elementary (parton-parton) ampli- 
tude instead of Martin's real-part formula and the selec- 
tion of more suitable parametrizations for the free pa- 
rameters. With this, we achieved a good (but limited) 



description of the main physical observables in the accel- 
erator energy domain. Based on these satisfactory results 
we extrapolate our parametrizations to cosmic-ray ener- 
gies in order to investigate total cross sections. 

In this report we present in some detail all the under- 
lying aspects of the work, related to both the technical 
matters and to the general ideas and physical interpreta- 
tions. 

The material is organized as follows. In Sec. |l| we 
briefly recall the experimental status (accelerator and 
cosmic-ray data) on pp scattering and outline a survey 
of some theoretical results that show the present rele- 
vance of the multiple diffraction formalism. In Sec. Ill 
we describe the model, the improvements introduced, the 
predictions for the physical observables in the accelera- 
tor energy region and the extrapolations to cosmic-ray 
energies. Discussions, physical interpretations and criti- 
cal remarks concerning the main results are the content 
of Sec. ffvl. 



II. 



EXPERIMENTAL AND THEORETICAL 
CONTEXTS 



Elastic pp scattering is essentially characterized by 
three physical observables: differential cross section, total 
cross section and the p-parameter (ratio of the forward 
real to imaginary part of the scattering amplitude) . From 
these, others quantities may be derived such as the slope 
parameter and integrated elastic and inelastic cross sec- 
tions ^. 

Experimental data on these observables obtained and 
compiled nearly fifteen years ago, still remain the only 
source of information at high-energies extending up to 
^/s = 62.5 GeV, the CERN Intersecting Storage Ring 
(ISR) energy region [^jj^. Experimental information on 
pp total cross sections in the range of energy ^/s : 6 — 40 
TeV exists from cosmic-ray data on Extensive Air Show- 
ers However, since the proton cross section is ex- 
tracted from the proton-air cross section through phe- 
nomenological models the results are model depen- 
dent. Also, the analysis performed by the Akeno Collab- 
oration 0] was recently criticized by N.N. Nikolaev who 
claimed that the Akeno results have been underestimated 
by about 30 mb |^ and this is in agreement with the re- 
sults early obtained by Gaisser, Sukhatme and Yodh 
This is a central point in our work and we will discuss 
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this discrepancy in Sees. II] and IV. Anyway, before the 
CERN's LHC proton-proton coUider, cosmic-ray results 
remain the only source of information on pp interactions 
at the highest energies. 

From the theoretical point of view it is well known that 
a pure QCD description of the elastic hadron scattering 
has not yet been obtained. Perturbative QCD cannot 
be extended to the soft region and pure nonperturbative 
QCD is not able to predict scattering states. Although 
the bulk of experimental data has been successfully de- 
scribed by different models in different contexts |^], a 
widely established and accepted approach is still miss- 
ing. The same can be said of a pure QCD understanding 
of the subject. 

As stated before, we are interested in the interface 
phenomenology/ nonperturbative QCD. For this we will 
limit this introductory discussion to some ideas directly 
related to this subject. 

Despite the above difficulties, there has been in recent 
years important progresses in the framework of nonper- 
turbative QCD, concerning soft processes. Looking for 
a microscopic approach able to explain phenomenolog- 
ical results, Landshoff and Nachtmann associated the 
Pomeron to the exchange of two abelian gluons with 
modified propagators ||]. Following this analysis, Nacth- 
mann extended the approach to the case of non-abelian 
gluons in a functional representation of scattering ma- 
trix elements and eikonal approximation With this, 
the quark-quark amplitudes are associated with a gluonic 
correlation function. However, since physical observables 
are connected with hadron-hadron amplitudes, a central 
problem concerns the construction of these amplitudes 
from the elementary ones. 

One way to treat this problem is based on the stochas- 
tic vacuum model The hadronic amplitude is con- 
structed through scattering amplitudes for Wilson loops 
in Minkowski space-time (loop-loop scattering), leading 
to gauge invariant amplitudes |[Tl|] . With this nonpertur- 
bative framework general characteristics concerning to- 
tal cross sections and slope parameters in high-energy 
hadron-hadron scattering may be described [ |ll| . 

Another way to construct the hadronic amplitude from 
the elementary one is by means of Glauber's multiple 
diffraction theory and this is the point we are inter- 
ested in. The approach is based on the impact parameter 
and eikonal formalisms as follows: Assuming azimuthal 
symmetry in the collision of two hadrons A and B, the 
impact parameter formalism connects the eikonal x ^^id 
the elastic hadronic scattering amplitude by [T^ 

F{q, s)=i bdb[l - e'''^''''^]Jo{qb) = i<l- eM''^'^ >, 
Jo 

(1) 

where = —t is the four-momentum transfer squared, b 
the impact parameter, Jq the zero-order Bessel function 
and the angular brackets denotes the symmetrical two- 
dimensional Fourier transform. In the first order Glauber 



multiple diffraction theory the eikonal is expressed as the 
Fourier transform of the products of the hadronic form 
factors, Ga and Gb, by the averaged elementary (parton- 
parton) amplitude /, |]l3| , namely: 



x{b,s)=<GAGBf> 



(2) 



This approximation means that at any time one con- 
stituent of a hadron interacts with only one constituent 
of the other hadron. This corresponds to the generalized 
form of the Chou-Yang model once a well-defined 
Fourier transform for the elementary amplitude is as- 
sumed ]l5| . 

The importance of this phenomenological framework 
in the search for connections between experimental data 
and calculational schemes in nonperturbative QCD has 
been recently expressed in the work by Grandel and 
Weise . Making use of the Dosch and Kramer ansatz 
for the gluonic correlation function [0 the authors cal- 
culated the elementary (parton-parton) amplitude f{q) 
in the high-energy approximation. Then, introducing 
a monopole parametrization for the unknown hadronic 
form factors in (||) the hadronic amplitude (0) could be 
calculated, leading to the differential cross section. With 
this, the authors obtained a satisfactory description of 
the differential cross section data for pp and pp elastic 
scattering at ISR and CERN Super Proton Synchroton 
(SppS) energies, in the region of small transfer momen- 
tum (g2 < 2.0 Gel/2) ^ 

Despite all these developments we see that the con- 
nections between nonperturbative QCD approaches and 
the bulk of physical observables still depends strongly on 
phenomenological models and ad hoc parametrizations. 
In this sense we understand that the role of a constant 
feedback on phenomenological information is crucial for 
present and future developments. 

In particular we see that the multiple diffraction for- 
malism may be a powerful tool in the test of suitable 
parametrizations. Simultaneously it may contribute with 
the search of reliable calculational schemes at differents 
levels of the theory (e.g. the stochastic vacuum model). 
We stress also the importance of the energy dependences 
which may be extracted from parametrizations for the 
form factors and elementary amplitudes. Since the the- 
oretical approaches (nonperturbative) mentioned above 
make use of asymptotic energy limits, the full increase 
of the total cross section, the shrinkage of the diffraction 
peak, the energy dependence of the p-parameter and even 
possible residual differences between pp and pp scatter- 
ing at the highest energies have not yet been explained. 

Based on all these considerations, in the next sections 
we will investigate elastic pp scattering in the context of 
a multiple diffraction model. 



III. MULTIPLE DIFFRACTION MODEL 
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A. Previous approach and Martin's formula 

From a phenomenological point of view, multiple 
diffraction models are currently identified by each par- 
ticular choice of parametrizations for the form factors, 
GajGb and elementary amplitude, / in Eq. (||) [ p^ . 
With this, the hadronic amplitude (|l|) may be investi- 
gated and then, in principle, the physical observables re- 
ferred to in the beginning of Sec. || 



atot{s) = ATiIm{F{q = 0,s)}, 



Re{F{q^Q,s)} 
Im{F{q = 0,s)}' 



(3) 



(4) 



(5) 



In a previous approach, elastic pp and pp scattering 
were investigated in a multiple diffraction model through 
the following parametrizations for the form factors and 
elementary amplitude [|l7|,|l8[ 



G = (l + ^)-^(l + ^) 



q 



f^iC- 



1 



1 



(6) 



(7) 



where a^,/3^,a^ and C are free parameters. We intro- 
duced here a small change in the previous notation of 
the elementary amplitude which will be suitable in what 
follows. 

The justifications for the above parametrizations were 
extensivelly discussed and explained in and and 
will not be reproduced here. We only recall that Eq. (Q) 
is in agreement with the results from model-independent 
analysis for the dynamical part of the eikonal as obtained 
by Buenerd, Furget, and Valin and Carvalho and 
Menon ||2l|. Discussion on this subject may be found in 
and llJ. 

With the parametrizations (^) and the eikonal in 
(2) is purelly imaginary and so the hadronic amplitude 
( 1 ) . In the previous approach the real part of this ampli- 
tude was estimated through Martin's formula |23 



ReF{q,s) = —[pq^ImF{q,. 



(8) 



using the experimental p value at each energy. With this 
approach a satisfactory description of experimental data 
on differential and integrated cross sections was obtained 
0,[9|j24|. 

However, a crucial point concerns the use of the above 
formula and this deserves some discussion. First, as 
derived by Martin the formula holds only for values of 



the momentum transfer infinitesimally small and at the 
asymptotically high-energy regime [ P3| . The formula 
may also be derived through the geometrical scaling hy- 
pothesis and in this case its applicability should be 
limited to the ISR energy region. Corrections to Mar- 
tin's formula were introduced by Hcnzi and Valin |2^ ] 
and numerical analysis from fits of experimental data by 
Kundrat and Lokajicek puts serious limits in its appli- 
cability concerning momentum transfer [ pTf . This result 
however has been recently criticized by Kawasaki, Mae- 
hara and Yonezawa who present results favouring the 
applicability of the formula in the whole region of the 
momentum transfer with the data available [ p8| . 

Despite this controversial aspect, a serious problem 
with the use of Martin's prescription is that p{s) is an 
input parameter at each energy and so cannot be inves- 
tigated. This led us to try a different procedure in the 
determination of the real part of the hadronic amplitude 
and in the next section we introduce a possible solution 
for the problem. 



B. Complex elementary amplitude 

In the context of the multiple diffraction formalism, 
Eqs. (0) and (||), associated with a complex hadronic 
amplitude, Fab , one should expect a complex elementary 
(parton-parton) amplitude 

f{q, s) = Ref{q, s) + ilmf{q, s). 

In this sense, the approach discussed in the last section 
corresponds to the assumption of Ref{q, s) — and 



Imfiq,s)^C- 



1- % 
1 + 4' 



(9) 



Lacking both theoretical and experimental information 
about the elementary phase we will assume, as a first 
approximation, a proportionality relation between real 
and imaginary parts at each energy p9| 



Refiq,s) = X{s)Imf{q,s), 



(10) 



where A(s) is a free parameter. With this assumption the 
eikonal in (||) may be expressed by 



x{b,s) = {x + t)nib, s) 

where, for the proton-proton case, 

r!(6,s)-<G2/m/(g,s)>, 



(11) 



(12) 



with G given by and Imf{q, s) by (^). With this, the 
real and imaginary parts of the hadronic amplitude read 



Re{F{q, s)} =< e^^^^''''*) sin(Af^(6, s)) >, (13) 



Im{F{q, s)} =< 1 - e-^'-''^'^ cos(Af](6, s)) > . (14) 



3 



Substituting parametrizations (||) and (||) for G and 
Imf in ( |l2|) , the "opacity" 17(5, s) is analytically evalu- 
ated. Then, Eqs. (|l|) and lead to the differential 
cross section (||). 

As a first test of the ansatz (10) we calculate the 



imaginary part of the hadronic amplitude, using yet the 
parametrizations from the previous approach jr7| , p^ ; 
then the real part was evaluated through both, Martin's 
prescription (||) and the ansatz (|l^). For pp elastic scat- 
tering at y/s — 52.8 GeV (largest interval in momentum 
transfer with data available, [^) the differential cross sec- 
tion is well described for A = 0.055 and we display both 
results in Fig. ^ 

We see that the predictions for the real part are similar 
in both cases: they present two zeros (change of sign) 
and its contributions to the differential cross sections are 
important only in the dip region. 

Since the imaginary part of the amplitude presents a 
zero and the real part of the Martin formula is obtained 
by the derivative (|^) , the contribution of this part is out 
of phase and then the differential cross sections does not 
vanish in the dip region With the assumption ( p^ ) 
for the elementary amplitude the same effect is obtained 
in the hadronic amplitude due to sin and cos terms in 
( p^ ) , (|lj) , and the result is in qualitative agreement with 
the predictions of zeros (change of sign) from dispersion 
relation ]3C| ]. 

Based on this satisfactory result with the complex el- 
ementary amplitude, we investigated the simultaneous 
descriptions of all physical observables referred before, 
Eqs. (||), (Q) and (||). To our knowledge this has never 
been achieved through geometrical or multiple diffraction 
models and we will return to this point later. 

Trying global descriptions for pp elastic scattering, be- 
sides the use of the complex elementary amplitude, we re- 
analysed the fits and parametrizations, improving some 
aspects of the previous approach. We discuss in detail 
the central points in what follows. 



C. Fits of experimental data 

As explained in the last two sections, our approach 
has only five free parameters: two associated with the 
form factor, and /3^, and three with the elementary 
amplitude, C, and A. 

We analysed 7 sets of pp experimental data above 10 
GeV (Table and the fits were performed only of the 
differential cross section data and p parameter data 

at each energy. The fit procedure consists of two 

steps: 



1. Taking A = in Eg. ( |11| ) the hadronic ampli- 
tude, Eqs. ( |l3|) and (|lj), is purely imaginary. For 
this case we determined the values of the parame- 
ters C, , and that reproduce the differential 
cross section data at each energy so as to present 
the zero at the dip position. 



2. With the values of the above four parameters as 
input we then calculated the value of A that repro- 
duces the experimental p value at each energy. 

Step 1 comes essentially from the previous approach and 
has been discussed and explained in detail in |l^ . 

With the above procedure a satisfactory description of 
p and da/dq^ experimental data was obtained with two 
constant free parameters 



= 8.20 GeV^, 



= 1.80 Gey^ 



(15) 



and only three parameters depending on the energy: 
C{s),a^{s) and A(s). The values are shown in Table || 
for each set analysed. 

In order to obtain a formalism able to make predic- 
tions to other energies we then proceed to investigate 
parametrizations for the data displayed in Table ^ 



D. Parametrizations as a function of the energy 

The choice of suitable and consistent parametrizations 
is a crucial point, mainly concerning extrapolations. We 
will discuss this aspect in detail in what follows and 
also in Sec. IV. We first analyze the dependences of C 
and q;~^, which determines the imaginary part of the 
hadronic amplitude (and so the total cross section) and 
after the dependence of A (associated only with the real 
part of the amplitude) . 

The values of C and a^^ from Table | are displayed 
in Fig. 1^ and we see that both increase with the energy, 
presenting positive curvatures. 

Experimentally, total cross sections grow like [Ins]", 
n = 1 or 2, at and above ISR and there is indication of 
"qualitative" saturation of the Martin-Froissart bound, 
n = 2.2 ±0.3 1^,0. Also, from gauge field theory, lowest 
order cross sections for particle production (unitarity) 
present In s terms . Based on these facts and from 

the behaviour shown in Fig. ^ we introduced fits through 
polynomials 



JV 



On [In 



n=0 



So 



which is different from our early parametrizations 
[p^-p^. Both sets of points are statistically consistent 
with polynomials of second degree [^6| and through lin- 
ear regression we obtained [^ 

C{s) = 14.3 - 1.65[ln(s)] + 0.159[ln(s)]2, {GeV'^) 

(16) 

^7 = 2.57 - 0.217[ln(s)] + 0.0243[ln(s)]2, [GeV''^) 

(17) 
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and we took So = 1 GeV^. In Fig. | we show the above 
parametrizations together with the fit values from Table 
H Physical interpretations and critical remarks concern- 
ing param etrizations (|l^ ) and ( [l7| ) will be presented in 
Sec. fVBl. 



E. Model predictions and experimental data 



In the case of parameter A the choice of a suitable 
parametrization demands further discussions. Empiri- 
cal analysis of the influence of A in the hadronic am- 
plitude showed that its behaviour is similar to that of 
the p parameter. That is, if A increases (decreases) also 
p increases (decreases) and A = at the same energy 
value where p = 0. We will return to this point in Sec. 
IV C. Since there is no experimental information on p{s) 
above 62.5 GeV we have, in principle, serious limitations 
in the choice of parametrizations. That is, in order to 
make extrapolations to high energies in our strictly phe- 
nomenological approach, we should infer the p behaviour 
or investigate limiting cases. With this last possibility in 
mind, we recall the general belief that, above ^/s ^ 100 
GeV, p{s) has a maximum and then goes asymptotically 
to zero through positive values 0|. However, how fast 
this happens depends on model assumptions. In order 
to test these possibilities we considered two different ap- 
proaches to zero, which could be understood as some kind 
of limiting cases, that is, slow and fast convergences. The 
point as we shall show is to investigate the influence of 
these assumptions in the description of the experimental 
data. 

In Fig. ^ we display the values of A from Table |. Based 
on this behaviour and on the above considerations we 
introduce the following general parametrization: 



X{s) = 



Ailn{s/sn) 



1 + A2[ln{s/so)\ + A:i[ln{s/sQ)Y ' 



(18) 



In this formula so controls the point where A(s) (and 
p{s)) reaches zero and Ai, z = 1,2,3 the maximum and 
asymptotic behaviour. As the two limiting cases we took: 



Case I: Ai= 6.95 x 10"^^ A2 = 0.118, 
A3 = 1.50 X 10-2 

Case 2 : Ai = 9.08 x 10-^, A2 = 0.318, 
A3 = 1.70 X 10-1" 



(19) 



(20) 



and So = 400 CeV'^ in both cases. Figure || shows these 
parametrizations up to ^/s = 10^ CeV . 

With the parametrizations for C{s),a^'^{s) and A(s) 
, Eqs. (jl^), (p^ ) and (18), respectively (cases 1 and 2), 
and result (|15|) for the remaining parameters 0^ and a^ , 
all free parameters are complet ely de ter mined. Through 
the formalism described in Sees. [II A and III B| , Eqs. (6), 
(3) and (O) to (p^), the three physical observables (3), 
(I) and ^ may be predicted. 



As explained in Sec. HI C our flts were performed only 



on differential cross section and p parameter data in the 
interval yfs = 13.8 — 62.5 GeV. In this session we first 
check the predictions in the interval ^/s = 10 — 100 GeV 
(accelerator energy region) and then display the extrap- 
olations up to 10^ GeV (cosmic-ray energies and future 
accelerators). 

• Accelerator energy region 

Figures ^ and || show the model predictions (cases 1 
and 2) for the pp differential cross section, p parameter 
and total cross sections, together with the experimental 
data available. We observe that in this interval, cases 1 
and 2 are distinguishable only for p{s) and ij^^^{s) and be- 
low -^s ~ 15 GeV. With the exception of the diffraction 
minimum (dip) at the highest ISR energies, the agree- 
ment with the experimental data is quite good. In Sec. 
[V C we will return to this overestimation of the differ- 



ential cross section in the dip region. 
• Extrapolations to higher energies 

Assuming that our parametrizations hold at higher en- 
ergies we calculate the predictions for p(s), a^^^{s) up to 
y/s = 10^ TeV (cosmic-ray region) and for the differential 
cross sections at 10, 15 and 20 TeV (future LHC). 

We show in Fig. ^ the predictions for p(s) in cases 1 
and 2 (Fig. ^ and the experimental data available. We 
understand that case 1 should be the most reliable from 
a conservative point of view. The similarities between 
p(s) and A(s) (Fig. ^) will be discussed in Sec. IV Cj . 

The predictions for the total cross section arc shown in 
Fig. ^, together with accelerator data (Fig. ||) and results 
from cosmic-ray experiments which, due to the existing 
discrepancies, we briefly review in what follows. 

The information available on proton-proton total cross 
section, crfo^, from cosmic-ray air showers are obtained 
through the p-air inelastic cross section, cr^j^j^-^. How- 
ever, either the determination of the (Tpfaiy or the rela- 
tion between a^fl^.i^ and crfoj are model dependent |^ . In 
the detailed analysis of data from Fly's eyes experiment, 
Gaisser, Sukhatme and Yodh (GSY hereafter), estimated 
the limit 

crfo^t > 130 m6 ai ^/s - 30 TeV. 

Making use of the Chou-Yang relation between a^^^ and 
the slope parameter, they calculated Q 

afP = 175t^? mh at ^^ 40 TeV. 

More recently, based on analysis of the extensive air 
shower, the Akeno Collaboration presented results in the 
interval ^/s : 6 — 25 TeV. In particular they found [Q 

all^ = 133 ± 10 m6 ai ^^ = 40 TeV, 
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which is incompatible with the GSY result. However, 
Nikolaev showed that (J^j^^^j. inferred by the Akeno group 
should be identified with an absorption cross section and 
that this originates an increase of ^ 30 mb in the Akeno 
results for af^^ . From the Nikolaev analysis |^ , 

crf^t = 160 - 170 mbat yG= 40 TeV, 

which is in agreement with the GSY calculation. 

Although the Akeno results are usually referred to in 
the literature, the analysis either by Nikolaev or by GSY 
seems correct and both are incompatible with Akeno. 

All these informations are displayed in Fig. |^, together 
with accelerator data and our predictions in cases 1 and 
2. We see that the model predictions are in complete 
agreement with Nikolaev and GSY results and this im- 
plies in a faster increase of the total cross section than 
usually believed. In particular we predict 



and 



a^J^t = 147 mbatyG= 16 TeV, 



af^t = 176 mb at y/^ = 40 TeV. 



Further discussion on the subject will be presented in 



Sec. [VA 



Concerning the differential cross section, we calcu- 
late the predictions in the region to be reached by the 
CERN's LHC, which are displayed in Fig. |. The results 
present no dip but a shoulder at \t\ ^ 0.5 GeV^ and at 
\t\ > 4 GeV^ the curves decrease smoothly with no other 
structures. 



IV. DISCUSSION 

The model described presents three parameters that 
depend on energy: a^,C and A. In this last section we 
discuss some critical points and physical interpretation 
concerning the parametrizations and predictions. 



A. High-energy extrapolations and total cross 
sections 

Our parametrizations ( ]l6| ) and ( p7[ ) for C(s) and 
q;~^(s), respectively, were based on experimental infor- 
mation only in the interval 13.8 < ^/s < 62.5 GeV and 
in both cases the points indicate an increase with the en- 
ergy (Table ||, Fig. ||) . A crucial question concerns this 
limited region since different explicit parametrizations, 
statistically consistent with the set of points, may differ 
arbitrarily when extrapolated to higher energies. 

However, constraints on the choice of parametrizations 
may be found through physical information available and 
this played an important role in our approach. As men- 
tioned in Sec. HID, our strategy was based on the gen- 
eral assumption of the expected In s - behaviour of soft 



processes and on the reasonable hypothesis of polyno- 
mial functions on Ins. With these constraints, we have a 
linear system in the free parameters of the polynomials 
(with the exception of the assumed value sq = 1 GeV) 
and the statistical solution is unique leading to the 
forms ( [l6| ) and (|l^) which describe the points quite well. 

On extrapolating the predictions for the total cross 
section to cosmic-ray energies, our results agree very 
well with both analyses by Nikolaev |^ and by Gaisser, 
Sukhatme and Yodh § (Fig. 0). Since the approaches 
by these authors are totally independent of the consider- 
ations and assumptions we have made, and, as far as we 
know, there is no criticism concerning their results, the 
agreement shown in Fig. ^ suggests a real increase of the 
pp total cross section faster than generally expected. In 
particular we predict 

crf^t = 91.6 mb at ^=1.8 TeV, 

which is higher than even the CDF result for pp, a^^^ = 
80.03 ± 2.24 mb |||. This seems to favour the "Odderon 
hypothesis" |Q , a problem which still "remains entirely 
open both from the theoretical as well as from the exper- 
imental point of view" |El| . 



B. Blackening and expansion 

The parameter coming from the hadronic form fac- 
tor (P) is associated with the radius defined by 



i?2(s) = -6 



dG{q, s) 



and from (||) 



R{s) = (0.483) [ 



1-] 



(21) 



Through parametrization ( |l7|) for a ^ and 0^ 



1.80 GeV"^, Eq. 
as shown in Fig. 



151) , the radius increases with the energy 
yTWe can then interpret the parameter 
expansion effect" 



The parameter G corresponds to the "absorption con- 
stant" in the Chou-Yang picture 42 1 and is associated 
with the number of constituent partons in the context of 
the Glauber approach ||l^,|l^ . It then controls the "black- 
ening effect" coming from the absorption. Our results, as 
in the previous approach |l^,|l^, mean that hadrons be- 
come blacker and simultaneously larger as the energy in- 
creases, in agreement with the "BEL behaviour" (Black, 
Edgie and Large) found by Henzi and Valin [^ . 

These effects may be investigated through the be- 
haviour of the Inelastic Overlap Function, Gin{b,s), 
which is calculated from the unitarity condition in the 
impact parameter space |44] 



2i?er(5,s) = |r(6,s)|2 + G,„(6,s), 



(22) 
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where T(b, s) is the Profile Function, the Fourier trans- 
form of the hadronic amplitude (|l|): 

r(6, s) = 1 - e*x(''''*^ (23) 

In Fig. |l^ we show Gin{b, s) as function of the energy 
and for some fixed values of the impact parameter b. In 
the region 10 to 100 GeV (part (a)) we observe the si- 
multaneous increase of Gin at all values of the impact 
parameter. In part (b) we show the central {b = 0) and 
peripheral {b = 1/m) regions extrapolated to ^/s = 10^ 
GeV. We observe that in the central region the curvature 
becomes negative above y/s = 100 ~ 200 GeV and in the 
peripheral region, above -y/s ~ 10'^ GeV. The black disc 
limit {Gin = 1) however seems very far to be reached, 
i.e., much higher than 10^ GeV. 

Since C(s) and a^(s) control the blackening and ex- 
pansion effects, respectively, the dimensionless quantity 
Ca^ gives information on the influence of each effect 
as function of the energy. Figure |l^ shows the predic- 
tions obtained by means of parametrizations (|l^) and 
( p^ ) up to 10^ GeV, together with the values from fit 
(Table ||). We observe a minimum at ^/s ^ 30 GeV, a 
change of sign in the curvature (becomes negative) above 
y/s ~ 10^ GeV and the asymptotic limit value ^ 6.5 
(from the parametrizations). 



C. Forward real-to-imaginary ratios of the partonic 
and hadronic amplitudes 



the predictions overestimate the differential cross section 
in the dip region at the highest ISR energies, leading to 
a limited description of the set of physical observables. 
With the exception of this point, all the predictions are 
in agreement with the experimental data (Figs. ^and|^). 
Since the total cross section is calculated from the imag- 
inary part of the hadronic amplitude our novel results at 
cosmic-ray energies are independent of the effect in the 
dip region. 

We observe that the ansatz (|l^) may be formally equiv- 
alent to some other geometrical and multiple diffrac- 
tion models characterized by complex eikonals p0| , ^5| -^ . 
The novel point here was to treat this hypothesis explic- 
itly and investigate its consequence in the context of the 
Glauber approach, as, for example, the strong correlation 
between p{s) and A(s) seen in Figs. |^ and ^. 

To our knowledge, simultaneous and complete descrip- 
tions of cross sections (differential and total) and the 
p parameter, still remain an open problem in geomet- 
rical and multiple diffraction models. For this reason, 
we hope, our results may bring insights for further and 
deeper developments. 
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An essential characteristic of the Glauber multiple 
diffraction formalism is to connect elastic scattering cross 
sections for composite particles (originally nuclei and af- 
ter nucleons) with the scattering amplitudes of their indi- 
vidual components (originally nucleons and after partons, 
respectively) ||l^ , [l3| , ^ . In this context, the assumption 
( p^ ) of proportionality between real and imaginary parts 
of the elementary (parton-parton) amplitude means that, 
in particular, 



A(.) 



Refiq = 0,s) 
Imf{q^ 0, s)' 



i.e., at the partonic level A(s) plays the same role as p{s) 
at the hadronic level, Eq. (||). The similarities referred 
to in Sec. |IIID| between A(s) and p{s) may be seen by 



comparisons of Figs. ^ and^ (cases 1 and 2). 

The hypothesis that A does not depend on the momen- 
tum transfer is a very simple one and has been used here 
only as an ansatz. Despite this, from Fig. ^ the resulting 
contribution of the real part of the hadronic amplitude to 
the differential cross section is in qualitative agreement 
with the results obtained through Martin's prescription. 

On the other hand, in the context of our approach, 
the limitations of this simple ansatz appear when we try 
simultaneous descriptions of cross sections (differential 
and total) and the p parameter. As we show in Fig. 0, 
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FIG. 1. Results for the elastic pp differential cross section 
at ^/s = 52.8 GeV. The imaginary part of the arnplitude 
was calculated using previous parametrizations J|l^,[l8| and 
the real part through (a) Martin prescription (W^ and (b) 
the ansatz (|lo|), with A = 0.055. In both cases it is shown 
the contributions from the real part (dotted), imaginary part 
(dashed) and complex amplitude (solid). Experimental data 
on dajdt are from p| and the p value in Martin's formula 



FIG. 2. Values of the free parameters C and a ^, from Ta- 
ble |, that fit experimental data (circles) and parametrizations 
through ^ and (|l|) (solid line). 

FIG. 3. Values of the free parameter A, from Table ^ that 
fit experimental data (circles) and parametrizations through 
([|) in the cases (|l^ (dot) and (|^ (dash). 

FIG. 4. Model predictions for the differential cross sections 
in cases 1, Eq. (p^, and 2, Eq. ( |2c| ) (indistinguishable) and 
experimental data ( |3l) for ^/s = 13.8 and 19.4 GeV and ^ 
for the other energies). Curves and data were multiplied by 
factors of 10*^. 
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FIG. 5. Model predictions for the p parameter (left) and 
total cross section (right) in cases 1, Eq. (^^ (solid), and 2, 
Eq. (po|) (dash) and experimental data ( |^ and |Q for p and 
and for the total cross section). 

FIG. 6. Predictions for the p parameter in cases 1 and 2 
and experimental data (Fig. |^. 

FIG. 7. Predictions for the proton-proton total cross sec- 
tion and experimental informations: accelerator data pUs^ 
(crosses), Akeno (circles), Nikolaev (triangles), GSY 
limit at 30 TeV || (T) and GSY result at 40 TeV [| (square). 



FIG. 8. Predictions for the pp differential cross section at 
10 (solid), 15 (dash) and 20 TeV (dot). 

FIG. 9. Radius calculated through Eq. (|l]) with 
= 1.80 GeV^ and parametrization (|l^ for a"^. 



FIG. 10. Inelastic Overlap Function calculated from the 
eikonal through Eqs. (||) and 

FIG. 11. Predictions for the dimensionless quantity Co? 
(see text) and points from fit (Table Q). 

TABLE I. Values of the free parameters from fits of pp 
diferential cross section and p data at each energy. 



C(s) 

[GeV) (Gel/-^) 
13.8 9.970 

19.4 10.050 

23.5 10.250 
30.7 10.370 

44.7 10.890 

52.8 11.150 
62.5 11.500 



a-2(s) 


A(s) 




2.092 


-0.094 


4.77 


2.128 


0.024 


4.72 


2.174 


0.025 


4.71 


2.222 


0.053 


4.67 


2.299 


0.079 


4.74 


2.370 


0.099 


4.70 


2.439 


0.121 


4.72 
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